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If the lightest supersymmetric particle (LSP) is Higgsino-like, the thermal relic density is lower 
than the observed dark matter content for a LSP mass in the sub-TeV region. We outline constraints 
arising from the Fermi Gamma-Ray Telescope data and LSP production from gravitino decay that 
must be satisfied by a successful non-thermal Higgsino scenario. We show that in a generic class 
of models where anomaly and modulus mediated contributions to supersymmetry breaking are of 
comparable size, Higgsino arises as the only viable sub-TeV dark matter candidate if gravitinos 
are heavy enough to decay before the onset of big bang nucleosynthesis (BBN). The correct relic 
density can be obtained via modulus decay in these models. As an explicit example, we consider a 
modulus sector in effective field theory (D = 4, N = 1 supergravitiy arising from type IIB KKLT 
compactification) . Within this class of mirage mediation models, heaviness of the gravitino forces 
a sub-TeV Higgsino LSP and gives a Higgs mass around 125 GeV. In this example, the constraints 
from direct detection experiments are also satisfied. 



I. INTRODUCTION 

Supersymmetry not only stabilizes the Higgs mass 
against quantum corrections, it also provides a candi- 
date for dark matter. In i?-parity conserving models the 
lightest supersymmetric particle (LSP) is stable, hence 
a dark matter candidate. The lightest neutralino, which 
is a mixture of Bino, Wino and Higgsinos, is the most 
suitable dark matter candidate with the prospect for de- 
tection in various direct and indirect searches. 

In this work we point out that a comprehensive solu- 
tion to the cosmological gravitino problem motivates the 
dark matter to be Higgsino-like. Gravitinos heavier than 
O(40) TeV have a lifetime shorter than 0.1 s and de- 
cay before the onset of big bang nucleosynthesis (BBN). 
This results in a considerable relaxation as the gravitino 
abundance will not be subject to tight BBN bounds Q]. 

In effective supergravity, the masses of the Bino and 
Wino are sensitive to the mass of the gravitino m 3 / 2 0, 
and in particular, for m 3 / 2 > 40 TeV, one typically has 
Bino and Wino masses above TeV in type IIB modulus 
mediation models. On the other hand, the Higgsino mass 
depends on the fi parameter, which can be reduced by 
anomaly mediated contribution to supersymmetry break- 
ing. As a result, if we demand that the dark matter 
particle has a mass in the sub-TeV region, the Higgsino 
becomes a more natural candidate. 

If the lightest neutralino is pre-dominantly Higgsino, 
with mass in the sub- TeV region, the annihilation rate 
is typically larger than the nominal value (cr ann t>) = 
3 x 10~ 26 cm 3 s _1 , thus resulting in an insufficient ther- 
mal relic abundance t 3J. A natural way to obtain the cor- 
rect dark matter relic density is to consider non-thermal 
sources of Higgsino production. 

We consider scenarios where Higgsino dark matter is 
non-thermally produced by a late decaying modulus [I]. 



We find that for annihilation rate to be compatible with 
bounds from the Fermi Gamma-Ray Telescope [7], the 
modulus decay should reheat the universe to a temper- 
ature Td ~ O(GeV). An additional requirement is that 
the branching ratio for modulus decay to the gravitino is 
<; 0(1O~ 5 ), so that the decay of gravitinos thus produced 
does not lead to dark matter overproduction. 

As an example of the modulus sector, we consider the 
standard scenario of KKLT compactification [5] , with the 
Kahler modulus reheating the universe around 1 GeV. 
Within this framework, for appropriate values of the rel- 
ative contributions of anomaly and modulus mediated 
contributions, Higgsino emerges as the dark matter can- 
didate. The annihilation rate is consistent with the Fermi 
bounds, and the correct relic density is obtained by non- 
thermal production. The Higgs mass m/j ~ 125 GeV 1 is 
also satisfied in this scenario, and we find that it actually 
requires the gravitino mass to be in the cosmologically 
safe region. Moreover, the spin independent scattering 
cross section is consistent with the latest bounds from 
direct detection experiments [8]. 

Within this specific example, however, decay of the 
gravitinos that are directly produced from modulus decay 
overproduces dark matter. This is a direct consequence of 
the couplings between the modulus and the hclicity ±1/2 
components of the gravitino, which are in turn set by the 
underlying Kahler geometry of the effective D — 4, N — 
1 supergravity theory. We summarize a set of geometric 
conditions in the modulus sector that are sufficient to 
ensure consistent non-thermal Higgsino dark matter as 
outlined above. 



1 Recent experiments at the LHC have provided strong hints of a 
Higgs- like particle at ~ 125 GeV [6]. 
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We note that apart from the purely cosmological mo- 
tivations shown in this study, the Higgsino also emerges 
as the LSP within the framework of Natural Supersym- 
metry as discussed in jS], [TO], [3]. 

The paper is organized as follows. In Section [XXJ we re- 
late the cosmological gravitino problem with the prefer- 
ence for Higgsino dark matter. In Section [lll[ we outline 
the conditions that must be satisfied by any successful 
scenario of non-thermal Higgsino dark matter. In Section 
IV, we work out an explicit example of a non-thermal sce- 
nario. In Section [Vj we outline the general constraints 
on an effective modulus sector in order to avoid overpro- 
duction of gravitinos. We conclude the paper in Section 
VI. 



B. Higgsino LSP 

To obtain a TeV scale spectrum in the observable sec- 
tor for such heavy gravitinos, one needs to consider mod- 
els where there is a hierarchy between the gravitino and 
the other supcrpartncr masses. We consider modulus me- 
diation, where supersymmetry is broken by a gravitation- 
ally coupled modulus in the hidden sector. An example 
of a broad class of models that provide the required hi- 
erarchy is in compactihcations of Type IIB string theory, 
with the following ingredients: (a) a Kahler modulus T a 
stabilized by non-perturbative effects (6) complex struc- 
ture moduli stabilized by fluxes (c) visible sector on D7 
branes [12] ■ In such cases, the gaugino masses obey 



II. COSMOLOGICAL GRAVITINO PROBLEM 
AND HIGGSINO DARK MATTER 



In this section, we discuss the cosmological gravitino 
problem, and argue that requiring dark matter in the 
sub- TeV range makes Higgsino a natural candidate due 
to an interplay between modulus and anomaly mediation 
contributions. 



A. Cosmological Gravitino Constraint 

The decay width of a particle <fi, which may be a 
modulus or the gravitino, with couplings of gravitational 
strength to the visible sector fields is 



m 3/2 
ln(M P /m 3 / 2 ) 



(3) 



When anomaly mediated contributions are subdomi- 
nant, the Bino is the LSP. However, for m 3 / 2 > 40 TeV, 
the Bino mass is typically above 0(1) TeV as seen from 
Eq. ([3]) . We will demonstrate this in our explicit example 
later. 

For the Higgsinos, there is an additional freedom. 
Starting with equal modulus mediated contributions, 
anomaly mediation lowers the mass of the gluino, while 
increasing the mass of the Bino and Wino. In this case, 
one has 



M 3 : M 2 : M x 
~ (l-0.3a) 53 2 



(1 



0.la)g 2 2 



(1 + 0.66a) ff ? , (4) 



(1) 



where c depends on the couplings of the decaying field. 
For moduli fields, we typically have c ~ 0.1 — 1. For grav- 
itinos, c can be computed explicitly since supersymmetry 
fixes the couplings of the gravitino to the visible sector. 
One has a maximal value of c ~ 1.5 in this case [TTj . 

The decay occurs when H ~ Tj,, with H being the 
Hubble expansion rate of the universe. Modulus decay 
reheats the universe to the following temperature 



T d ~ (5 MeV) c 



1/2 



/ 10.75 

V 9* 



100 TeV 



3/2 



(2) 



where g* is the total number of relativistic degrees of 
freedom at T d (g t = 10.75 for T d > 0(MeV)). 

Gravitinos that have a lifetime shorter than 0.1 s de- 
cay before the onset of BBN and avoid any conflict with 
its successful predictions. Such alifetime corresponds to 
T d > 3 MeV, which requires that m 3/2 £ O(40) TeV 
from Eq. 



where Mq is the modulus mediated contribution at the 
GUT scale, a = m 3 / 2 /A/o ln(Mp/m 3 / 2 ) denotes the rel- 
ative strength of anomaly and modulus mediated con- 
tributions, and <7;l,2,3 are the gauge coupling constants. 
[T3l The limit a — > corresponds to vanishing 

anomaly mediated contribution and Bino-like LSP. 

The value of the Higgsino mass parameter /1 depends 
on the low-scale value of m 2 H , which is mainly driven 
by the gluino mass j!5j . Increasing a lowers the gluino 
mass, which in turn lowers the low-scale value of m? H 
due to the top Yukawa coupling. On the other hand, the 
Bino and the Wino become heavier. 

Thus, if we demand the dark matter candidate mass 
to be less than 0(1) TeV, Higgsino LSP is preferred by 
the cosmological gravitino problem, provided that the 
anomaly mediated contribution to the soft masses com- 
petes with the modulus mediated contribution. 

The mass hierarchy between the scalar masses and the 
gravitino mass is more model dependent, and depends 
on the curvature properties of the underlying Kahler 
manifold. Stop masses of 0(1) TeV, preferred by a 125 
GeV Higgs mass in the MSSM, are obtained in models 
with heavy gravitinos where the suppression is similar to 
Eq. §}. 
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III. SCENARIOS OF NON-THERMAL 
HIGGSINOS PRODUCTION 

In this section, we consider non-thermal production of 
Higgsino dark matter. Standard relic density calculations 
predict thermal underproduction of Higgsinos for masses 
less than about a TeV [3]. Thus, one is motivated to 
study non-thermal scenarios where enhancement of the 
relic density occurs naturally. Non-thermal production 
is inevitable if the modulus <fi that participates in super- 
symmetry breaking, or any other modulus, reheats the 
universe at a scale below the dark matter freeze-out tem- 
perature Tf ~ m x /25. 

The dark matter relic density from modulus decay is 
given by 



5 x 10" 



1 GeV\ 3 x 10~ 26 cm 3 s" 1 /T £ 



where (cr ann v)f is the annihilation rate at the time of 
freeze-out. 

The Higssino DM mainly annihilates into heavy Higgs 
bosons, W bosons and t quarks via S'-wave annihilation 
if m x has necessary phase space for these particles to be 
produced. The S'-wave nature of the annihilation implies 
that annihilation rate at the freeze-out time is essentially 
the same as that at the present time. The latter is con- 
strained by the gamma ray flux from dwarf spheriodal 
galaxies [7J: 

(ffannw) / <, 10~ 25 cm 3 s" 1 m x = 100 GeV , 

K„nv) / < 3 x 10~ 24 cm 3 s" 1 m x = 1 TeV. (6) 

In order to obtain the correct DM abundance, see 
Eq. jjj), one therefore needs to have: 



T d £ 0.4 - 1.6 GeV 



100 GeV - 1 TeV . (7) 



For Td ~ O(GeV), the corresponding modulus mass is 
found from Eq. ^ to be 



0(1000) TeV, 



(8) 



with the exact value depending on the decay modes of 
the modulus. 

Another non-thermal source of Higgsinos production is 
gravitino decay. Since gravitino decay occurs at a tem- 
perature -C O(GeV), and dark matter annihilation rate 
must satisfy the Fermi bounds (6j), annihilation is very 
inefficient at this time. As a result, the density of Higgsi- 
nos produced from gravitino decay is therefore the same 
as that of the gravitinos. Therefore we must have: 



^^5x 10- 10 



1 GeV 



(9) 



Gravitinos are produced via thermal and non-thermal 
processes in the early universe. Modulus decay dilutes 
gravitinos that were produced in the prior epochs (e.g., 



during inflationary reheating) by a huge factor. Thermal 
gravitino production after modulus decay is highly sup- 
pressed due to the low decay temperature Td ~ O(GeV). 

However, gravitinos can also be produced directly from 
modulus decay 4> —> GG. The density of gravitinos thus 
produced is given by (^3/2/s) = Br 3 / 2 (3Xd/4m0), where 
Br 3 / 2 is the branching ratio for <f> GG process. From 
Eqs. ([2]), we then find 



"3/2 



5 x 10" 



100 TeV 



1/2 



Br 



3/2 



(10) 



For the typical value of given in (JS|) and 100 GeV < 
m x < 1 TeV, Eqs. ( 9|10 ) yield the following absolute 
upper bound: 



B1-3/2 < 10" 



(11) 



Any successful scenario for non-thermal Higgsino produc- 
tion from modulus decay must satisfy this limit. 



IV. NON-THERMAL HIGGSINO DARK 
MATTER: AN EXAMPLE 

As an explicit example, we take the case of mirage 
mediation in the context of KKLT compactification [5]. 
Working in D = 4, N = 1 effective supergravity, the 
superpotential of the modulus sector consists of a flux 
term that fixes complex structure moduli, and a non- 
perturbative piece that fixes the Kahler modulus. The 
Kahler potential is given by K = -3 In (T + T) + (T + 
T)~ nm <I><I>t 5 where $ denotes matter fields and n m are 
the modular weights. The input parameters fixing the 
GUT scale masses are m 3 / 2 , ct, n m , and tan/3. 

For our case study, we choose n m = 1/2 for all matter 
fields and tan/3 = 50. The general conclusions hold for 
other values of n m and tan/3. 

The scalar spectrum has the suppression given in 
Eq. Q with respect to the gravitino mass. It is instruc- 
tive to note that when the stops are themselves hierar- 
chically related to the gravitino as in this case, ~ 125 
GeV is compatible with heavy gravitinos that decay be- 
fore the onset of BBN. Since the one-loop correction to 
the Higgs mass depends logarithmically on m 3 /2, a little 
heavier Higgs is preferred by the cosmologically safe re- 
gion. We plot the dependence of the Higgs mass on the 
gravitino mass in Figure [T] 

For given n m = 1/2 and tan/3 = 50, scanning over 
0.1 < a < 1.6, and to 3 / 2 > 40 TeV, one finds that for 
LSP mass below ~ 1 TeV, the LSP is always a Higgsino. 
We plot fj, against LSP mass in Figure [2] This result 
is independent of the choice of tan/3. Since the gaugino 
mass does not depend upon n m , this conclusion is also 
independent of n m . 

We also plot the spin independent scattering cross sec- 
tion crjQ—p for various values of the gravitino mass in 
Figure [3j Since larger gravitino mass is correlated to a 
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larger heavy Higgs (H) mass in this model, cr^o_ p 
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FIG. 1: Higgs mass versus gravitino mass. We choose n m = 
1/2 and tan/3 = 50, and scan over a and m 3 /2- Heavy graviti- 
nos decaying before the onset of BBN are typically compatible 
with a Higgs mass above 120 GeV. A similar behavior is ob- 
tained for other values of tan/3 and n m 
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comes smaller as m 3 / 2 increases. Again, this is compat- 
ible with cosmologically safe region and mt ~ 125 GeV. 
The current bound on the cross section is ~ 2 x 10~ 9 pb 
for a dark matter mass of 55 GeV [S], and relaxes as dark 
matter mass increases. 

Table [i] depicts a few benchmark points of consistent 
non-thermal Higgsino dark matter in this model. In the 
table, we show the annihilation rates at present (<7 ann t;}o 
and at the time of freeze-out (cr ann v)f. We see that for 
all of these points (<J ann v)o satisfies the Fermi bounds. 

Dark matter annihilation at the freeze-out occurs 
mostly through S'-channel and a coannhilation compo- 
nent. The latter arises due to the fact that masses of 
the second lightest neutralino and chargino are close to 
LSP mass m^i . The dark matter content of the universe 

AO 

is obtained by multiplying the (a ann v) / by Tf/T^. We 
have taken c = 0.4 to calculate T^, which is the leading 
order value appearing in the decay width of the modulus 
in this particular example, in the limit of a = 1 corre- 
sponding to zero dilaton-modulus mixing in the gauge 
kinetic function. But c can also be ~ 1 depending on 
relative contributions of the modulus and dilaton in the 
gauge kinetic function, and Td can be 0(1 — 2) of its 
central value. With this taken into account, it is seen 
that Tf/Td is in the right range to yield the correct dark 
matter content of the universe. 

We also show the value of cr^o_ p for these points, which 
are well allowed by the experimental data. The table also 
includes masses of the gluino and stops, and the Higgs 
mass rrih for these points. 
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FIG. 2: fi versus LSP mass. For sub-TeV LSP, the dark 
matter is always a Higgsino. 
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FIG. 3: Spin independent scattering cross section versus grav- 
itino mass. For values of TO3/2 that are compatible with 125 
GeV Higgs, see Fig. 1, cr^i_ p saT, i snes the experimental data. 



V. GENERAL CONDITIONS FOR 
SUPPRESSING GRAVITINO PRODUCTION 

In the KKLT model discussed above, the partial de- 
cay rate for -> GG is T 3/2 = m^/2887rM|. Then, 
after using Eq. Q, we find that Br 3/2 - 0(lfr 2 ). This 
implies that gravitino decay will overproduce Higgsinos 
by 3 orders of magnitude in this model, see Eq. (10 1. 



The main reason for obtaining such a large Br 3 / 2 is that 
modulus decay to helicity ±1/2 gravitinos is not helicity 
suppressed in the KKLT model [T6] . 

The problem can be overcome if the modulus (j) does 
not dominate the energy density of the universe when it 
decays. In such a case, the right-hand side of Eqs. (10) 
and (111 will be multiplied by and , respectively, 
where/0 is the ratio of the energy density in <f> to the 
total energy density of the universe at the time of de- 
cay. For /<£ < 10 -3 , the abundance of gravitinos will be 
suppressed to safe levels. 

Alternatively, one can seek conditions for suppressing 
gravitino production from modulus decay (for example, 
see [H]). Here we briefly outline general conditions for 
such a suppression, and stress that the main ingredients 
of a successful scenario for non-thermal Higgsino produc- 
tion presented above should also hold in cases where the 
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gravitino production is suppressed. 



TABLE I: Some benchmark points of non-thermal Higgsino dark matter for mirage mediation model in the context of KKLT 
compactification. The input parameters are a and m 3/ / 2 - The modular weights are fixed to be n m = 1/2, and tan/? = 50. All 
masses are in GeV. 
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m 3 / 2 


mh 




X9. 


m s 




m t 2 


(fannW)o (cm 3 /s) 


(cr a nnW)/ (cm 3 /s) 


T f /T d 


o^o_ p (pb) 


1.49 


143 • 10 3 


123.5 


248.4 


250.8 


3828 


2441 


2781 


1.49 ■ 10" 2b 


1.63 • 10"^ 


~ 6 


5-10" lu 


1.46 


200 • 10 3 


124.5 


258.9 


260.6 


5536 


3564 


3991 


1.38 • 10 -25 


1.52 • 10~ 25 


~ 3.4 


1.4- 10 -10 


1.44 


232 • 10 3 


125 


306 


308 


6505 


4197 


4677 


1.01 ■ 10~ 25 


1.01 • 10~ 25 


~ 3.2 


8.9- HT 11 



The decay of a modulus to other fields depends on the 
interaction terms in the Lagrangian, and the requirement 
for suppressing decay to gravitinos will be reduced to 
a set of constraints in the effective theory. To have a 
more concrete demonstration of what kinds of constraints 
may emerge, we choose to work in effective supergravity, 
with a modulus coupling to the visible sector through the 
gauge kinetic function. This is the scenario in the class 
of Type IIB models discussed above. 

In general, one can consider a scenario with multiple 
moduli <pi, with the decaying modulus appearing in the 
gauge kinetic function. The normalized eigenstates <p n 
are given by 



(<f>)i 



E 

3 



Gij {l~ n )j 



(12) 



where the Cy are eigenvectors of the matrix K^ 1 d 2 V. 
For simplicity, we will assume diagonal Cij with entries 
Ci. The partial widths for modulus decay to gauge fields, 
gauginos and helicity ±1/2 gravitinos are 



No 



1287T 



\ 2 m <t>i 



288tt 



mi /2 M 2 



(13) 



where G = K + log \W\ 2 is the Kahler function. 

Under suitable choices of the Kahler potential, the re- 
quired condition Br 3/ / 2 ~ 10~ 5 may be obtained. Simi- 
larly, the decay temperature of the modulus may be ob- 
tained in terms of the Kahler potential and superpoten- 
tial. We refer to [T7] for more details. 

For a single modulus, the branching ratios to gauge 
bosons and gauginos are roughly equal, and the branch- 
ing to the gravitino needs to be suppressed, leading to 
the condition 



WI3/2 



\GJK~j 



1/2 



10" 



(14) 



For the KKLT example, the above quantity is 0(1), 
which leads to overproduction of gravitinos. However, 



in a more general scenario, one can suppress this ratio 
to the required levels by a suitable choice of and 
vacuum expectation value of <f>. This does not necessar- 
ily affect the existence of other conditions for success- 
ful non-thermal Higgsino production, such as compara- 
ble anomaly mediated contributions, or a modulus in the 
correct mass range. Moreover, the scalar masses depend 
on the holomorphic bisectional curvature of the plane (in 
tangent space) spanned by the scalars and the super- 
symmetry breaking modulus, and this is not necessarily 
changed by a shift in the metric Kj,z. One can therefore 
expect to have a viable non-thermal scenario with the 
Higgs mass mh ~ 125 GeV, while suppressing gravitino 
production from modulus decay. We leave the detailed 
exploration of these issues for future work. 



VI. CONCLUSIONS 

Considering dark matter in the sub-TeV range, ther- 
mal freeze-out underproduces Higgsino-like LSP. It is 
possible to enhance the relic density using non-thermal 
mechanisms of dark matter production. The enhance- 
ment, however, needs to obey the constraints from the 
gamma ray flux from dwarf spheroidal galaxies and the 
dark matter content of the universe. Moreover, there 
should be no overproduction of dark matter at any later 
stage, for example by the decay of gravitinos. 

In this paper, we demonstrated these ideas in a generic 
class of mdoels where anomaly and modulus mediated 
contributions to supersymmetry breaking are compara- 
ble. Interestingly, we found that within this class of mod- 
els, heavy gravitinos that are not subject to BBN bounds 
force the Higgsino as the only viable dark matter candi- 
date in the sub-TeV range. We considered an explicit 
example of mirage mediation model in D = 4, N = 1 
supergravitiy arising from type IIB KKLT compactifica- 
tion, where the modulus decay provides the non-thermal 
origin of Higgsino-like dark matter. The large gravitino 
mass is helpful to yield rat around 125 GeV in this model 
and satisfy the constraints arising from the dark matter 
direct detection experiments. We also discussed the gen- 
eral conditions to avoid the overproduction of LSP from 
gravitino decay in such scenarios. 
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